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ABSTRACT 



Felix Bloch’s 1928 article made a prediction concerning the dynamical behavior 
of electrons in a solid, subject to a uniform, static electric field. This aspect of his work, 
as later clarified by Zener, showed that electrons accelerated by an electric field in a 
periodic potential, under the right conditions, would oscillate. A theoretical debate as to 
the existence of this phenomenon has been ongoing since Bloch’s proposal. One of the 
most controversial consequences of this prediction is that an electron undergoing Bloch 
oscillations would radiate. The controversy on the theoretical analysis was due to the 
great difficulty in systematically and reliably treating interband transitions by analytical 
methods based on the time-dependent Schrodinger equation for independent electrons. In 
this thesis, we numerically solve the time-dependent Schrodinger equation to show that 
electrons accelerated by in an electric field in periodic structures do undergo Bloch 
oscillations and other dynamic behavior. By accurately modeling this phenomenon, we 
hope to gain a better understanding of it in hopes of using it in future applications as a 
stable source of Terahertz (THz) radiation. 



VI 



TABLE OF CONTENTS 



I. INTRODUCTION 1 

II. THEORETICAL WORK 3 

A. FELIX BLOCH 3 

B. CLARENCE ZENER 4 

C. W.V. HOUSTON 6 

D. GREGORY WANNIER 8 

E. J. ZAK 9 

F. A. RABINOVITCH 10 

G. J.N. CHURCHILL AND F.E. HOLMSTROM 11 

H. J.B. KRIEGER AND G.J. IAFRATE 12 

I. L. ESAKI AND R. TSU 14 

J. D. EMIN AND C.F. HART 14 

K. M. LUBAN AND A. BOUCHARD 15 

III. NUMERICAL METHODS 17 

A. TIME DEPENDENT SCHRODINGER EQUATION (TDSE) 17 

1. Discretization of the Wave Function 18 

2. Discretization of the Hamiltonian 18 

B. SOLUTIONS OF THE TIME DEPENDENT SCHRODINGER 

EQUATION 21 

1. Power Series Expansion 21 

2. Magnus Expansion 22 

3. Explicit Algorithm for Solving the TDSE 23 

4. Cayley Hamilton Approximation 24 

C. SOLVING THE CAYLEY APPROXIMATION 25 

1. Solution by Inversion 26 

2. L-U Decomposition 25 

vii 



IV. RESULTS AND CONCLUSIONS 29 

A. NON LOCALIZED WAVE PACKET IN GENERAL POTENTIAL 30 

B. LOCALIZED ELECTRON IN GENERAL POTENTIAL 32 

C. ELECTRON REFLECTIONS AT BOUNDARIES 33 

D. KRONNIG-PENNY SIMULATIONS 34 

1. Non localized Electron in a Superlattice 34 

2. Effect of Varying Repeat Distance 35 

3. Effect of Varying Barrier Thickness 36 

4. Effect of Varying Electron Mass 37 

E. SUMMARY 38 

APPENDIX: CAYLEY METHOD IN ANSI C 39 

LIST OF REFERENCES 47 

INITIAL DISTRIBUTION LIST 49 



viii 



I. INTRODUCTION 



In 1928, Felix Bloch [Ref. 1] wrote a groundbreaking work discussing the 
quantum mechanical nature of electrons in solids. In this seminal work, Bloch 
demonstrated what is widely known as “Bloch’s theorem,” which establishes the 
mathematical form of the electron wave function of a crystalline solid. We will discuss 
Bloch’s theorem later in this thesis. Bloch’s theorem is the basis for our understanding of 
the all-important energy bands in solids. In addition to establishing Bloch’s theorem, a 
less well known result from Bloch’s 1928 article is a prediction concerning the dynamical 
behavior of electrons in a solid, subject to a uniform, static electric field. This aspect of 
his work, as later expanded by Zener [Ref. 2], showed that electrons accelerated by an 
electric field in a periodic potential, under the right conditions would oscillate. The 
frequency of oscillation is given by the formula a> B = eFa/ft , where eis the charge of an 
electron, F is the electric field amplitude, a is the lattice repeat distance, and ft is 
Planck’s constant h divided by 2 ji . This phenomenon has come to be known as Bloch 
oscillations. A theoretical debate as to the existence of this phenomenon has been 
ongoing since Bloch’s proposal. One of the most controversial consequences of this 
prediction is that an electron undergoing Bloch oscillations would radiate. The 
controversy on the theoretical analysis was due to the great difficulty in systematically 
and reliably treating interband transitions by analytical methods based on the time- 
dependent Schrodinger equation for independent electrons. 

With the advent and refinement of semiconductor technology, it has been finally 
possible to test the predictions of Bloch’s proposal. In 1970, Esaki and Tsu [Ref. 3] 
proposed using a semiconductor superlatrice, which would provide the needed periodic 
structure on a large enough scale in order to conduct tests to generate and detect Bloch 
oscillations. According to Esaki and Tsu [Ref. 3], the oscillating electron in a 
semiconductor superlattice would provide a source of Terahertz (THz) radiation. It was 
only in 1992, however, that Bloch oscillations and the associated Terahertz radiation were 
conclusively detected in semiconductor superlattices [Ref. 4]. These recent four wave 



mixing experiments confirm the presence of photons emitted as a result of the oscillating 
electron, but noted the oscillations died between 1 and 15 oscillations [Ref. 4]. More 
recently, Luban and co-workers[Ref. 5] modeled the effect of an electron in a 
semiconductor superlattice using computer techniques, and showed Bloch oscillations 
exist at Terahertz frequencies, but did not show the decay as seen in real world 
experiments. The advantage of a purely numerical approach to solving the time- 
dependent Schrodinger equation is that one explicitly works with the full Hamiltonian, 
and not some approximation as in analytical solutions, thus including all interband 
transitions in the simulation. 

This paper will cover the theoretical debates and arguments surrounding the issue. 
Also, it will develop the theory behind Bloch oscillations. Lastly, it will cover the 
numerical techniques involved in solving the Time Dependent Schrodinger Equation, and 
present the results of the simulations. 
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II. THEORETICAL WORK 



In this chapter, we summarize and discuss the works of major contributors on the 
subject of Bloch oscillations. 



A. FELIX BLOCH 



In 1928, Felix Bloch wrote a ground breaking paper [Ref. 1] which until recently 
was the subject of great debate. In it, he used quantum mechanics to model the effects of 
electrons in periodic (atomic) structures. In this landmark paper, he proved what is now 
known as Bloch’s theorem, which shows that the mathematical form of the eigenvector of 
the time independent Schrodinger equation must be: 

( 2<1 ) 

where 

“nA x ) = u nA x + a ) ( 2 - 2 ) 

and n labels the energy bands and is known as the band index. Also, the parameter A: is 
the allowed wavevector of the electron and is restricted to the Brillouin zone associated 
with the particular lattice structure of the crystal. The associated energy eigenvalue 
E n (k ) is called the band structure function. 

Bloch showed [Ref. 1] that the wave function of any particle in a periodic 
potential can be represented as a sum of orthogonal functions that have the periodicity of 
the lattice structure. Bloch [Ref. 1], starting with the time dependent Schrodinger 
equation with an external electric field applied: 



, 2m i2m 

V 2 ^(x, 0 - — (V(x) - eFx)y(x,t) + — ■ -3^ = 0 
n~ n dt 



(23) 
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